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C\ '. Abstract 
C3 : 

We show how observations of the perturbation spectra produced during 
■ inflation may be used to constrain the parameters of general scalar-tensor 

theories of gravity, which include both an inflaton and dilaton field. An 
interesting feature of these models is the possibility that the curvature per- 
turbations on super-horizon scales may not be constant due to non-adiabatic 
perturbations of the two fields. Within a given model, the tilt and relative 
\Q • amplitude of the scalar and tensor perturbation spectra gives constraints on 

. the parameters of the gravity theory, which may be comparable with those 

Q\ , from primordial nucleosynthesis and post-Newtonian experiments. 

cr 

bp. 

The most convincing explanation for the flatness, isotropy and homogeneity of the ob- 
served universe is the inflationary scenario [Q. Moreover, the most compelling evidence for 
this model is the prediction of a nearly scale-invariant distribution of Gaussian perturbations. 
If these are indeed the origin of the perturbations observed in the microwave background 
sky, and of the initial inhomogeneities from which galaxies formed, then they could provide 
our earliest glimpse of the physics of the early universe and, in particular, of the effective 
theory of gravity at that time. 

In this paper, we consider the possible constraints that can be placed upon the al- 
lowed gravity theory during inflation. Precision tests of gravity in the solar system severely 
constrain the effective gravity theory today ||, while predictions from primordial nucleosyn- 
thesis have been used to restrict scalar-tensor deviations from general relativity since the 
radiation dominated era Our aim is to push back those limits to a still earlier epoch, 
i.e. inflation. 

We will do this within the context of general scalar-tensor gravity theories, which involve 
a massless dilaton or Brans-Dicke field These provide a well-defined class of theories 

against which to test the predictions of general relativity. Almost all attempts to produce 
a renormalizable quantum theory of gravity seem to include scalar fields non-minimally 
coupled to the spacetime curvature. The low-energy effective action of string theory, for 



I. INTRODUCTION 
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instance, involves a dilaton field coupled to the Ricci curvature ||. Scalar fields coupled 
directly to the curvature tensor appear in all dimensionally reduced gravity theories, and 
their influence on cosmological models was first seriously considered by Jordan [|7]]. Gravity 
Lagrangians including terms of higher order in the Ricci scalar can also be cast as scalar- 
tensor theories J8]|| with appropriate scalar potentials. 

Recently, Damour and Nordtvedt [jT0[] have pointed out that during a dust dominated 
cosmological era, scalar-tensor theories of gravity tend to approach general relativity at late 
times even in the absence of a potential for the Brans-Dicke field. This can be parametrized 
by the Brans-Dicke parameter uj, which determines the ratio of the scalar to tensor couplings 
to matter, and tends to infinity in the general relativistic limit. Damour and Nordtvedt 
realized that this would occur during most of the history of the universe, when the trace of 
the energy-momentum tensor drives the scalar field. However, this mechanism is not effective 
during the radiation dominated era |jll|] . Considering a wide class of theories where u is some 
arbitrary function of the gravitational coupling, Damour and Nordtvedt [JRJ calculate how 
far towards the general relativistic limit the universe might be expected to evolve (i.e. how 
large to becomes) after the radiation era. They find that it should be simply related to the 
number of e-foldings, N = ln(ao/a eq ) since matter-radiation equality. 

As most of the expansion of the universe to date occurred during the inflationary era, 
one expects the same approach to general relativity to occur during inflation |12|]. This has 
been recently discussed in the context of string theory by Damour and Vilenkin JL3|| . Thus, 
apart from setting the scene for the conventional hot big bang, by producing a spatially 
flat, isotropic, homogeneous (but slightly perturbed) metric, it is tempting to suggest that 
inflation may also produce general relativity as the low-energy effective theory of gravity, 



even if it started as a generic scalar-tensor theory at the Planck scale [14 



We show in this paper that this can indeed occur. We give the generalization to scalar- 
tensor gravity of the inflationary slow-roll parameters in general relativity and show that 
the vanishing of these parameters corresponds to the usual post-Newtonian limit of Einstein 
gravity. We calculate the scalar and tensor perturbations produced during inflation in a 
general scalar-tensor theory and evaluate the relative amplitude and slope of the spectra 
in terms of the slow-roll parameters. While it is difficult to obtain model-independent 
bounds on the scalar-tensor parameters, we will study a particular model of chaotic inflation 
for which strong bounds can be given that are comparable with the post-Newtonian and 
nucleosynthesis limits. 



II. SCALAR TENSOR GRAVITY THEORIES 



The scalar-tensor field equations are derived from the action 
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S = d xv— g 



f{(j>)R - -g ab <p, a <P, b - U{4>) + £ matter 



(2.1) 



where R is the usual Ricci curvature scalar and 16n f((j>) is the Brans-Dicke field. Thus the 
gravitational coupling strength (Newton's constant G in general relativity) is determined 
here by the dynamical variable /(</>). In the particular case of Brans-Dicke gravity, /(</>) = 
4> 2 /8uj, where uj is the Brans-Dicke parameter, and we recover general relativity in the limit 
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u) — > oo. More general scalar-tensor theories with different choices of f(4>) correspond to the 
case where w(f) is a function of /. 

The potential U(<p) is the generalization of the cosmological constant A in general rel- 
ativity. Perturbation spectra produced in models of inflation driven by a potential for the 



Brans-Dicke field have recently been discussed by Kaiser while the bounds on the al- 
lowed mass have been discussed by Damour and Vilenkin |13| and Steinhardt and Will [16|. 
Such a potential is often introduced to fix the value of Brans-Dicke field at late times, 
however we shall show that this is not necessary in order to attain general relativity as a 
cosmological attractor. In what follows we will leave f(4>) as a free function but consider 
only models in which U(4>) is zero. They correspond to flat directions in the scalar poten- 
tial of string effective theories |17[ . In the absence of non-perturbative effects, the dilaton 
remains massless. We will work in this approximation and will discuss possible extensions 
in future work. 

Matter is minimally coupled to the metric g ab and thus test particles follow geodesies in 
this frame, which we refer to as the Jordan frame. However it is often useful to write the 
action in terms of the conformally related Einstein metric 



g ab = e' 2a ^g ab , (2.2) 
where the conformal factor is defined as e~ 2a = 2k 2 f (0), in terms of which the action in 



Eq. ( p.l| ) takes the Einstein-Hilbert form with a fixed gravitational constant G = k 2 /87t, 
and the Brans-Dicke field can be written as a scalar field ip with a canonical kinetic term in 
the new matter lagrangian 

^matter = —^a^fi + ^ ^matter • (2.3) 

There are now explicit interactions between this field and the original matter fields whose 



energy-momenta are therefore not necessarily conserved with respect to g ab [IB,|I9|. We 
could define a dimensionless parameter [|l(J, a(ip) = n~ l (da / dip) , which characterizes the 
scalar-tensor theory: a 2 specifies the ratio of the dilaton and graviton couplings to matter. 
A given choice of function a(ip), or equivalently /(</>), determines a(ip). In particular, it is 
related to the Brans-Dicke parameter u by 

2a 2 = . (2.4) 

2^ + 3 v ; 

For a linear a(ip) = anip with constant a, we recover Brans-Dicke theory. [] 

Present day observational tests constrain the Post-Newtonian parameters 7 and (3 [0, 
written in terms of a and a' (if)) = n~ 1 (da/ 'dip) as 



1 Ouv notation coincides with that of Damour-Gibbons-Gundlach [20] and Starobinsky 
Yokoyama |21] for constant —2a = 7 = (5/2. 
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which are constrained by present-day solar system tests to be H (7 — 1| < 0.002 and 



\4(3 - 7 - 3| < 0.001 ||. Therefore, 

« 2 <5-10" 4 , 4a 2 |l + 2a'\ < 10~ 3 . (2.6) 

There are similar constraints from primordial nucleosynthesis [|J. These constrain possible 
variations of the Planck mass during and after the radiation dominated era. Our aim is to 
go beyond the radiation era and try to constrain possible deviations from general relativity 
during inflation. 

III. SCALAR-TENSOR INFLATION 

In this section we will analyze the classical evolution of the scalar fields during inflation. 
The inflaton field, a, minimally coupled in the Jordan frame, with a self-interaction potential 
V(a) gives an explicit matter lagrangian to consider in a scalar-tensor cosmology. The total 
matter Lagrangian in the Einstein frame including the Brans-Dicke field is then 

^matter = -"^V,^ " <f V^, b - 6^ V (a) . (3.1) 

We see by inspecting the potential term in the Lagrangian that o will evolve towards 
a minimum of V(a) while if) evolves towards a minimum of a(if>) or, equivalently, a zero of 
a(ifj). However, from Eq. (|2.4|) we see that a zero of a(if>) requires that uo — > oo. That is, 
general relativity will generically be an attractor during the cosmo logical evolution, if a(if)) 
possesses a minimum. 

The field equations for the fields a and xf> in a spatially flat Friedmann-Robertson- Walker 
metric are then 

a + 3H& = - 2ani)a - e 2a V\o) , (3.2) 
ip + 3Hip = - a/t(e 4a AV - e 2a a 2 ) , (3.3) 



K 
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H = -'-{^ 2 + e 2a cr 2 ) , (3.4) 



and the Hamiltonian constraint 
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F 2 = ^(V> 2 + e 2a <7 2 + e 4a 2y) . (3.5) 



The condition for inflation to occur in the Einstein frame \H \ < H is thus, see Eqs. ( |3.4| 



and ( |3.5| ), 

ijj 2 + e 2a & 2 <e 4a 4V(a) . (3.6) 
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A. Slow-roll inflation 



We will work in the slow-roll approximation in both scalar fields. In principle this is not 
a necessary constraint: one of the fields might roll quickly to the minimum of its potential 
and then the problem reduces to single field inflation, either the familiar chaotic inflation in 
general relativity (for ip = 0) or old extended inflation in Brans-Dicke (for & = 0). However, 



we would like to consider the more general case in which both fields slow-roll |23|-25|. In 
this case, the general field equations can be written as first-order equations, 

3H 2 ~ K 2 e 4a V(a) , (3.7) 

3Ha ~ - e 2a V\a) , (3.8) 

3HtP ~ - Aane^Via) . (3.9) 

Neglecting the other terms in the equations of motion amounts to the following assump- 
tions 

max{^ 2 , a 2 ^ 2 , e 2a a 2 } < e 4a V(a) , (3.10) 
\a\ < \H&\ and \$\ < \Hip\ . (3.11) 

Having written down first-order equations for the evolution of the fields we can turn the 
slow-roll assumptions based on values of the fields' derivatives into consistency equations in 
terms of the parameters of the theory: 

^4wT' ( 3 - 12 ) 

= 8a 2 (ip) , rfy = 4a'(^) - 16a 2 (if>) . (3.13) 

The consistency equations for slow-roll inflation are then {e~ 2a e a , e~ 2a \rj a \, e^, \r]^\} 1. 
The first two conditions are just the expected generalization to scalar-tensor gravity of 
the slow-roll conditions for an inflaton field in general relativity. Notice, however, that if 
a = Katp < during the subsequent evolution of the universe (i.e. a > 0), the conditions on 
e a and rj a are relaxed compared to the general relativistic case, where a(i[)) = throughout. 
The last two parameters are not present in general relativity and arise here by requiring 
that both the Brans-Dicke and the inflaton field slow-roll. 

We have defined the slow-roll parameters e and rj by extending the usual definition of 




these parameters for a single field [26[ to the separable potential for the two fields U(a } ip) = 
V(o)e 4a ^\ It is intriguing to note that the limit of vanishing slow-roll parameters for 
the ip field coincides with the general relativistic weak-field limit in the post-Newtonian 
parametrization of the scalar-tensor gravity theory ||, see Eq. ( |2.5| ). 

In calculating the rate of change of quantities with respect to different co-moving scales, 
it is useful to write down the relation between the number of e-foldings from the end of 
inflation and the values of the scalar fields, 

N = -f Hdt= K - [ -% = n 2 [ ^ da . (3.14) 



4 V a(ip) Ja e e 2a V'(a) 
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Our present horizon crossed outside the Hubble scale about 50-60 e-foldings before the end 
of inflation. In fact, the precise number depends logarithmically on the energy scale during 
inflation and the efficiency of reheating, and so is weakly model-dependent. 



IV. DENSITY PERTURBATIONS 

Inflation is the only known mechanism that solves the horizon and homogeneity prob- 
lems. However, the main observational constraint on inflationary models is the spectrum of 
density perturbations that they produce. Strictly speaking, observations of perturbations 
in the microwave background, or of the large-scale structure in our patch of the universe, 
only provides an upper limit on the level of density perturbations, which could perhaps be 
produced by some other source of inhomogeneities. Nonetheless, the apparently Gaussian 
and nearly scale-invariant nature of the perturbations are natural properties of perturbations 
due to quantum fluctuations of the inflaton field during inflation. 

In the case of a single slow-rolling field, only adiabatic perturbations are possible. Any 
fluctuation in the field must produce a fluctuation in the local curvature. However, in the 
presence of two coupled fields we must also consider the effect of isocurvature (or entropy) 
perturbations between the two fields. In particular this can lead to the breakdown of the 
usual assumption that the curvature perturbation is frozen-in on scales greater than the 
Hubble length (/c p h < H). It is important to allow for such effects if we hope to constrain 
the possible role of a variable gravitational coupling. 

Our calculations extend those of Starobinsky and Yokoyama |T| who considered the 



particular case of inflation in Brans-Dicke gravity. As we shall see, their results are readily 
extended to more general scalar-tensor theories. 

A. Perturbed field equations 

In this section we will consider the linear perturbations around the homogeneous back- 
ground fields, a {t) + 5cr(t, x % ), ip(t) + 8ip(t, x l ), with a perturbed metric 

ds 2 = - (l + 2$(t, ar*)) dt 2 + R 2 (t) (l - 2$(t, x 1 )) 5 ij dx i dx j . (4.1) 

We can study the evolution of each Fourier mode (whose physical wavenumbers we denote 
by k ph ) separately, since they decouple in the linear approximation. The perturbed field 
equations then yield the following expressions to first order 

8a + 3H5a + k 2 ph 5a + e 2a V"(a)5a = 4a<5> + -2e 2a V'{a)<S> 

- 2aK(&8ip + ^5<j) - 2a'(ip)K 2 8tp - 2an e 2a V'(a)6t/j , (4.2) 

£0 + 3H Sip + k 2 ph 5tfj + a'{tfj)n 2 (e 4a AV - e 2a a 2 ) 5^ + 2a 2 n 2 (e 4a 8V - e 2a a 2 ) 5tfj 

= 4^$ - 8a/te 4a V<$> - 2an (e 4a 2V\a)8a - e 2a a8a) , (4.3) 
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<£ + ah§ + (h + m 2 ) $ 



+ e 2a &5o - e 4a V'{a)5a + e 2a a 2 - e 4a 4F are ^ 



together with the energy and momentum constraints 



K 2 



#V> + e 2a o5o + e 4a V'(a)<^ + (e 2a a 2 + e 4a 4^) a« <^ 



(4.4) 



(4.5) 



(4.6) 



A very useful quantity for the study of perturbation spectra is the three-curvature of 



comoving hypersurfaces |27| , [28 



tt2 

if v 



(4.7) 



Combining Eqs. Q4.4|) — fl4.6|) , and using the equations of motion, we find an exact expression 
for the time variation of £, 



H 
H 



d e 2a a 2 -ip 



dt \e 2a a 2 + ?p 2 , 



+ C 




(4.8) 



where C = 2otKipe ia a i / (e 2a a 2 + ip 2 ) 2 is due to the frictional damping of the a field by if). 

In the single field case (where one of the fields is held fixed) the right-hand side of 
Eq. ( |4.8j ) vanishes in the limit /c p h — > 0, and thus ( remains constant on scales exceeding the 
Hubble length p8| . This allows one to determine the large-scale curvature perturbation at 
the end of inflation simply by equating it with the perturbation when the mode first crossed 
outside the Hubble scale. However this is true in general only for adiabatic perturbations 
and need no longer hold in the presence of two fields p9j^T| . 



This is due to the entropy perturbation pO 



t5S= H 



d 



'e 2a & 2 



^ 2 



dt \e 2a & 2 + ip 2 , 



C 




(4.9) 



The first term in the square brackets will be present whenever two fields are evolving but 
the second term, C, would not be present if both fields had standard kinetic terms. It is 
clear that for adiabatic modes 5ip/ip = bo jo (perturbations along the classical trajectory) 
( remains constant on super-Hubble scales, but in general the curvature perturbation at 
horizon crossing cannot be equated with that at the end of inflation, due to the non-adiabatic 
terms. 

Note that all the above perturbed equations are exact and we have not yet invoked the 
slow-roll approximation. In the following we shall solve for the evolution of ( to lowest order 
in the slow-roll parameters. 

We will assume that the curvature perturbation will be conserved on super-Hubble scales 
through re-heating and the subsequent radiation and matter dominated eras. Recently, it 
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has been emphasized |3T[] that the curvature perturbation ( is conserved at transitions in 
the equation of state across a boundary at a fixed energy density. This is the case when 
perturbations are adiabatic and the end of inflation must coincide with a particular energy 
density. However this is not necessarily the case in the presence of two fields. We can 
see from Eq. ( |4.8| ), that our assumption that only adiabatic perturbations exist at the end 
of inflation requires that the motion of one of the fields dominates. This will usually be 
the case, especially when inflation ends because one of the fields' kinetic energy becomes 
comparable to the potential energy. The slow-roll parameter e for one of the fields becomes 
of order unity while the other remains small. It is possible that both slow-roll parameters 
become large at the same point, but this seems to be unlikely a priori. 

The perturbation at the end of inflation can then be equated with that at re-entry if 
its subsequent evolution remains adiabatic. Any variation of the Brans-Dicke field after the 
end of inflation could invalidate this assumption. Nucleosynthesis limits @ suggest this does 
not occur at temperatures below about an MeV and in the absence of a potential for the 
Brans-Dicke field, ip = is a stable solution during a radiation dominated era for arbitrary 



scalar-tensor gravity theories 11]. However this could be altered by the presence of an 



explicit potential term and the consequences would require careful investigation. 



B. Short-wavelength limit 

For large values of fc p h 3> H we can neglect the potential terms in the perturbed field 
equations ( |4.3| ) and ( }4.2| ) and they reduce to those for massless fields (i.e., e~ 2a \r]^\, \r] a \ <C 1). 
Thus, to lowest order in the slow-roll parameters, the expectation values of the perturbations 
as they cross outside the Hubble radius (k p ^ ~ H*) are given by Gaussian random variables 
with (|5<r*| 2 ) = e~ 2a *H 2 /2k 3 , (\5ip*\ 2 ) = H 2 /2k 3 , where k is the comoving wavenumber. 
Note that, while the field ip is minimally coupled in the Einstein frame, the a field is 
minimally coupled in the Jordan frame and therefore it is the conformally transformed 
Hubble constant (to lowest order) that determines its amplitude at Hubble crossing JEJ. 

We shall denote the spectrum of a quantity A by Pa(^) = j0js (\A\ 2 ) , as defined 



33]. Thus we have 

2 



7V-<T 2 - (^j , (4.10) 
V st * (§ )' . (4,1) 



C. Long-wavelength limit 

For slowly varying (<3> <C HQ), long- wavelength (k^ — > 0) modes, to lowest order in the 
slow-roll parameters, the Eqs. (fl5l), (|4.2|) and (f4.3|) reduce to 



$~-2aK5xp-^^-^-8a, (4.12) 
3H5ip ~ Aa' {ip) k 2 e Aa V 5ip , (4.13) 



S 



3H6a ~ - e 2a (^y^J V 5a + 2a Ke 2a V'(a) 5^ . (4.14) 

Note that for constant a we recover Starobinsky and Yokoyama's results [21]. 

Using Eqs. (|3.9j ), the last two equations can be integrated to give the evolution of fluc- 
tuations in the scalar fields at long-wavelengths: 

4a 

5^ ~ Q x , (4.15) 

K 

5a^^ X ^{Q 2 -e- 2a Q 1 ) , (4.16) 

and thus 

$ ~ 8a 2 Q t - JL (Q 2 - e^Qi) = t^Q, - e a (Q 2 - e~ 2a Q x ) , (4.17) 

where Qi and Q 2 are constants of integration, chosen to coincide with those introduced 
by Starobinsky and Yokoyama pl| . It will be convenient to define a new constant Q3 = 
Qie~ 2a * — Q 2 , so that Qi and Q3 are independent Gaussian random variables whose values, 
for a given Fourier mode, are determined by the amplitude of 5a* and Sip* at horizon- crossing 
(when k ph = H*). Thus they have expectation values 

p 4a„ 4 y 

v - - • < 4 - 18 > 

„2a» 4 y 

v* = w ■ («») 

From Eqs. (|4.5p and (|4.7[ ) we have, during slow-roll [ 25| , pi |, in the long wavelength limit, 

ijj 2 + e 2a a 2 

As shown in Eq. (|4.8|) this expression need not be constant. Substituting in our results for 
the long-wavelength modes of the scalar fields, we have 



^ ^ + (e- 2a -e~ 2a *)e u )Q 1 + e a Q, 



If either of the scalar fields is fixed (e CT or e$ identically zero) then we recover the single field 
results where ( is constant (equal to Qi or Q 3 respectively). 

The spectrum of density perturbations at the end of inflation V^ e (k) can be computed 



from (4.21 



2a»\^e\ 2 / „e \2 



h 4+r j v ^[^i) <422) 



In Sect. [VI] we will study a particular model and give numerical results showing how and 
when the different terms dominate. 
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D. Gravitational wave perturbations 



In addition to the scalar curvature perturbations that give rise to density perturbations, 
tensor or gravitational wave perturbations can also be generated from quantum fluctuations 



during inflation [p^j . Since we have chosen to work in the Einstein conformal frame, we 
can use the standard results for the evolution of tensor perturbations of the metric. The 
two independent polarizations evolve like minimally coupled massless fields with a spec- 
trum 

V a = . (4.23) 

Graviational wave perturbations can contribute to the microwave background 
anisotropies only on the largest scales (scales larger than the Hubble scale at last-scattering, 
corresponding to about > 1° on the sky). Their contribution relative to scalar curvature 
perturbations is given by the ratio 

3 V 

R ~ . (4.24) 

The rapid decay of the gravitational wave anisotropies on smaller scales is their most dis- 
tinctive signature. In Sect. |VT| we will study a particular model and show how the ratio R 
changes with scale. 



V. OBSERVATIONAL CONSTRAINTS 

Having allowed for the possible evolution of the curvature perturbation ( on super-Hubble 
scales during inflation, we will now restrict our analysis to those cases where ( has become 
constant on observable scales by the end of inflation, i.e., entropy pertubations become 
negligible. This allows us to assume that ( remains fixed on super-Hubble scales until it 
re-enters the Hubble length at late times. We can then relate the curvature perturbation 
at the end of inflation to the density perturbation at re-entry during the matter dominated 
era, 

4 - 7^ ' (5 ' 1} 



following the notation of 33 



In any model of inflation, the amplitude of the density perturbations depends upon the 
magnitude of the potential energy density relative to the Planck scale, which is essentially 
a free parameter. We will concentrate upon the variation of the amplitude of the curvature 
perturbations with co-moving scale. At each point in the spectrum, the 'tilt' is given by 
the spectral index n s , where n s — 1 = d In 5 2 H jd In k . This can be evaluated within the 
slow-roll approximation, where the comoving wavenumber corresponds to a given scale at 
horizon crossing, dink ~ —dN*, and thus from Eq. (|3.14j) , 



n _ 1 _ d\n5 2 H _ _ /4c^^_ + e 2a *Vj d \ ^ g2 
dink \ k dip* k 2 V* da* J 
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A well-known result in general relativity, for slow-roll inflation with a single field, is that 



n s — 1 ~ — 6e* + 277*, a function solely of the slow-roll parameters at Hubble-crossing [35,33 
Because ( can evolve on super-Hubble scales during scalar-tensor inflation, ( e and thus b 2 u 
will in general depend upon parameters both at Hubble crossing and at the end of inflation. 
However, the latter will not change with co-moving scale. 

Large-angle microwave background experiments probe scales close to our observable hori- 
zon, which crossed outside the Hubble scale when N ~ 60 [Q. The COBE two-year data 



constrains the spectral tilt to be in the range 0.7 < n < 1.7 at the la level ||36|| , but future 
experiments should be able to constrain the tilt to within about 0.1 P7[ . However, from 
observations of galaxy clustering we might hope to recover the primordial density perturba- 
tions on scales down to about IMpc, which leave the Hubble scale at about iV ~ 50. This 
could provide a much more precise determination of the tilt, though it is currently limited 
by uncertainties in other cosmological parameters such as the value of the Hubble constant 
or the type of dark matter. While a tilt of 0.7 < n < 0.8 may have some advantages 
over n = 1 in an otherwise standard cold dark matter model, a tilt n < 0.6 is probably 
unacceptable p8|j33[| . 



The general expression for the tilt from Eq. ( |5.2|) is rather complicated. We will only 
attempt to evaluate it in two general regions of parameter space and then specialize it to the 
case of an inflaton with a generic chaotic inflationary potential, for which we give numerical 
results. 

Note that the tilt of the gravitational wave spectrum is just given by 

n g ~-l6al-2e 2a *C . (5.3) 

Unlike the approximate expressions for the scalar tilt which will be given below, this simple 
expression for n g is valid in the whole range of parameter space. Moreover, since both terms 
on the right-hand-side must be non-positive, one could in principle give a direct constraint on 
a 2 completely independent of the form of the inflaton potential. However, the measurement 
of this slope will be exceedingly difficult. Tensor perturbations do not contribute to structure 
formation and in many inflationary models the observable effect of gravitational waves is 



completely negligible [33 



The main constraint coming from gravitational waves will be their relative amplitude, 
given by Eq. (|4.24j) . If R becomes of order unity then, since independent Gaussian ran- 



dom variables add in quadrature, the amplitude of the scalar perturbations inferred from 
anisotropies of the microwave background on large scales is reduced by about 70%. As R 
increases, the allowed amplitude of scalar perturbations decreases, eventually becoming in- 
compatible with structure formation. It is the combined effect of a tilted spectrum and the 
gravitational wave contribution that proves such a strong constraint on models of extended 
inflation |35|[ . 

Finally, note that if n s < 1 then, on very large scales, the potential energy density 
relative to the Planck mass at Hubble crossing becomes so large that the amplitude of density 
perturbations is of order unity and the universe enters a self- reproducing regime, where the 
classical motion is dominated by quantum fluctuations P^,|S2]JT^]. Such inhomogeneities, 



even on super-horizon scales today, could be detected through the Grishchuk-Zel'dovich 
effect [^] unless they are on scales at least 500 times greater than our horizon |JD|. This 



represents another constraint on any model. Sufficient inflation in the classical regime thus 
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requires iV max > 66. 



A. Scalar-tensor extended inflation 

Let us first consider the case where the Brans-Dicke field evolution dominates that of the 
inflaton at the end of inflation, eS, 3> e%. Then, for scales crossing outside the Hubble scale 
near the end of inflation, we have ( e ~ Q\ and thus 5 2 H ~ (A/25)Vq 1 . This remains valid at 
scales for which (e^/e^) 2 ^> e~ 2a *e^/e*. It includes models of extended inflation ^j] where 
the field a is trapped in a metastable false vacuum so that e a = (for any r\ a > 0) and where 
( remains fixed on super-Hubble scales. But this result for 5h also includes perturbations 
for which e* > e^, where there will be significant evolution of the curvature perturbation ( 
when k ph < H during inflation. 

We then have 

n s - 1 ~ -16a 2 + 8c< -2e 2a *e* a . (5.4) 

When e* = this expression generalizes the well known result for extended inflation in 
Brans-Dicke models |£J to more general scalar-tensor theories. We see that, just as in 
general relativity, n s need not always be less than unity. For instance, we can produce a 
Harrison-Zel'dovich spectrum (n s = 1) by choosing a scalar-tensor theory where a' = 2a 2 
corresponding to a(ip) = — (1/2) ln(ip/ip e ). This is a particular realization of "intermediate 
inflation" 0. 



More generally, as e a is always non-negative, a lower bound on the tilt of the power 
spectrum then constrains the slow-roll parameters of the gravity theory, irrespective of the 
form of the inflaton potential. 

The relative contribution of tensor and scalar perturbations to the microwave background 
anisotropies is given by Eq. ( |4.24| ), which in this limit yields 



i?~10 2 a 2 . (5.5) 
An upper limit on this ratio then constrains a 2 independently of a'. 



B. Scalar-tensor chaotic inflation 

In the opposite limit, e e a 3> e%, in which the evolution of the inflaton dominates that of 
the Brans-Dicke field at the end of inflation, we find ( e ~ Q3 and thus 5 2 H ~ (4/25)Pq 3 . 
This result will hold for the last scales to leave the Hubble length during inflation. It remains 
valid on larger scales as long as 

-2a* 




l-e- 2 - + ^ (5.6) 



e 2a *e 



is satisfied, see Eq. (|4.22 ). In these limits, the curvature perturbation is always due to 
fluctuations in the inflaton field er, but there may still be evolution on super-Hubble scales 
due to the frictional damping by ip. We find from Eq. ( |4.21|) that Q — e 2a Q^i which coincides 



with the solution of ( ~ C(, given by Eq. ([4.8|) in this limit 
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Given the above result for 5^ we thus find that on sufficiently small scales the tilt will 
be given by 

Us _ 1 ~ e 2a, ( _ 6e * + _ ga 2 (5 7) 

Note that r\ a can be positive or negative and thus could lead to a positive spectral tilt |33| . [43| . 

The larger effective gravitational constant at early times (a* > 0) amplifies the tilt due 
to the changing shape of the inflaton potential, and its variation leads to an additional 
negative tilt. Any chance of constraining 0? from observations of the tilt is clearly limited 
by uncertainty in the form of the inflaton potential. In the simplest case of chaotic inflation 
driven by a polynomial potential, V(a) = Xa 2n /2n, there is a simple relation between e a , 7]^ 
and the value of a, 

( n \ 2n 2 . , 

This guarantees that —660- + 2r/ CT is negative and thus the slope of the density perturbation 
spectrum in Eq. ( |5.7|) is always n s < 1. In arbitrary scalar-tensor theories, a lower limit on 
the slope then places an upper bound on a 2 . 

The ratio between the scalar and tensor contributions to the microwave background 
anisotropies reduces to the usual general relativistic case . 

R~l2e 2a *e* a . (5.9) 

When the condition given in Eq. (|5.6| ) no longer holds and instead (1 — e~ 2a *) ^> e^/e*, 
we find that ( e — Qi- In this regime the results of Eqs. ( |5.4j ) and ( |5.5| ) apply. It is interesting 
to note that the naive calculation based on taking C — C* does in fact give the correct result 
and, even in our careful analysis, ( remains constant on super-Hubble scales. This is clearly 
seen in the Fig. 1, where we show the evolution of £ after Hubble crossing in the specific 
model of chaotic inflation discussed in the next section. 

This occurs despite the fact that at Hubble crossing the curvature perturbation is due 
to the ip field, ( ~ H5ip/ip, while by the end of inflation it appears as a perturbation in the 
a field, ( ~ H5a/a. This is a consequence of the coupling between the two fields and the 
dependence of 5a upon the evolution of 5ip seen in Eq. (|4.14j ). We do not expect this result 
to hold in general for two fields in general relativity. Moreover, for intermediate scales ( 
does evolve on super-Hubble scales. 



VI. NUMERICAL RESULTS 

In this section we will try to show the main features discussed above with a particular 
example that includes both regimes. We will choose the arbitrary origin of the field ip so that 
ifj e = at the end of inflation (when a e = 0) and assume that a(ip) can be approximated as 
a linear function during the latter stages of inflation (when observable scales cross outside 
the horizon), i.e., we take a Taylor expansion for a(ip) up to second order, as was done in 0. 

a(ip) = di nip + (Kip) 2 , (6.1) 
a(ip) = ai + a 2 Kip , (6.2) 
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where a 2 and a% are constants and ip is then given in terms of the number of e-foldings, 
using Eq. ( p.l4j ), as 



Kip 



Oi 
«2 



[exp (4a 2 N) - 1] 



Therefore we have 



and 



a 



a\ exp (4a 2 N) , 



-2a 



2k 2 /(0) ~ exp 



-i (exp(8a 2 iV) - 1) 
a-2 



(6.3) 



(6.4) 



(6.5) 



Note that at the end of inflation (N = 0) we have the present value of the gravitational 
coupling, and our parameters a\ and a 2 correspond to a and a' respectively, at the end of 
inflation. 

For a simple inflaton potential V(er) = Acr 4 /4, the slow-roll solution for a, satisfying 
e% = 1, is 



K 2 <7 2 ~ 



+ 



exp(a 2 /a 2 ) 
a 2 



E 1 (a 2 /a 2 ) - Ei(a\/a 2 ) 



(6.6) 



where E\(z) is the exponential integral function []44 |. 
It will be convenient to define some new variables 



x 



KG 



y = exp(-2a(V0) 



(6.7) 



whose classical evolution in the slow-roll approximation is shown in Fig. 2 for a 2 = 10~ 2 and 
a\ = 10~ 2 and a\ = 10~ 3 . Note that in both cases y, and thus the Planck mass, becomes 
essentially constant by the end of inflation. 

If we also introduce w = 1 — y + 8a 2 , then in the limit w 2 y*x 2 a 2 , where all starred 
quantities are to be evaluated at N — iV*, we satisfy the condition given in Eq. (|5.6| ) and 
we find Ce — Q^- This must hold for the last modes to leave the Hubble scale at the end of 
inflation, as y — > 1 and x — > 2 a/2, for 8a 2 -C 1. For larger scales, when w 2 ^> y*x 2 a 2 , we 
find ( e ~ Q\. This is clearly demonstrated in Figs. 3(a) and 3(b). 

The spectral tilt ( |5.2|) at a scale which crosses outside the Hubble length when N = N*, 
is given, at lowest order in the slow-roll parameters, by 



(2a 2 - a 2 + (2/y^xl)) w 2 + 2a 2 y,w, + a 2 (3 + y*x 2 a 2 ) 



wl + y*x 2 a 2 



(6.8) 



see Fig. 4. In Fig. 5 we have also computed the ratio of gravitational to scalar components, 

96a 2 (8a 2 + l) 2 



R 



wl + y*xlal 



(6.9) 



We thus recover the results of Eqs. ( p.7| ) and ( |5.9j ) for n s — 1 and R respectively on small 
scales, and on larger scales by Eqs. (|5.4j ) and ( |5.5| ). The variation of both the tilt and the 
ratio R with Hubble crossing epoch, iV*, is shown in Figures 4 and 5. 
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Observational constraints on R and the tilt of the perturbation spectrum can bound the 
values of the parameters a\ and ci2, thus constraining deviations from general relativity as 
far back as iV* e-folds before the end of inflation. Figure 6 is a contour plot showing n s for 
scales that left the horizon at iV* = 60, corresponding roughly to 6000 Mpc today, and thus 
the sort of scale constrained by observations of large-scale structure. We see that both a\ 
and ci2 must be very small in order for the tilt of the spectrum to remain close to the general 
relativistic value of n s ~ 0.95. This reflects the need to keep the Planck mass essentially 
constant to avoid large departures from the Harrison-Zel'dovich (n s = 1) spectrum. Figure 7 
shows similar results for the contribution of gravitational wave perturbations at the same 
scale, N* = 60. 



VII. CONCLUSIONS 

In this paper we have considered the constraints that may be placed upon the effective 
theory of gravity during a period of inflation in the early universe. We do this in the 
context of scalar-tensor theories, taking the coupling of the Brans-Dicke field to matter as 
an arbitrary function a(ip), and neglecting any explicit potential for the dilaton field. 

Present day observational limits on the variation of the Brans-Dicke field are expressed 
as bounds on the post-Newtonian parameters of the theory. We have shown that the general 
relativistic limit of these parameters coincides with the vanishing of the corresponding slow- 
roll parameters for the Brans-Dicke field during inflation. Slow-roll inflation already requires 
the scalar-tensor theory to be close to the general relativistic limit. The observed spectrum of 
density perturbations produced from quantum fluctuations in the inflaton and Brans-Dicke 
fields can then constrain just how large the deviation may be. 

A careful calculation of the curvature perturbation ( during inflation shows that some 
of the results, applicable to single-field inflation in general relativity, no longer apply. Due 
to the evolution of the two fields during inflation, there will be non-adiabatic perturbations, 
which can lead to the evolution of ( on scales larger than the Hubble length. Therefore, 
the amplitude of ( at re-entry can no longer be equated with that at the time the scale left 
the horizon. This is a general feature of inflation with two fields. However, in scalar-tensor 
theories there are regimes for which ( does remain constant outside the Hubble scale. In 
models where inflation ends as the inflaton field rolls to the minimum of its potential, we 
find two regimes for which ( at the end of inflation is equal to that at Hubble-crossing. It 
is only in an intermediate regime that the naive calculation breaks down. 

We give expressions for the spectral slopes, n s and n g , and relative amplitude, R, of 
the scalar and tensor perturbations produced, in terms of the slow-roll parameters of a 
general scalar-tensor theory of gravity. Observational bounds then place constraints on 
these parameters. A possible signature of scalar-tensor inflation is the breakdown of the 
consistency relations predicted in single-field inflation [EH . 



It is important to emphasize that our ability to make quantitative predictions relies on 
our knowledge of the inflaton potential. In our specific example of a chaotic inflation model, 
we find that if we constrain the slope of the spectrum to be n s > 0.6 we obtain bounds 
on a and a' at the end of inflation that are comparable with those from nucleosynthesis or 
solar-system tests. For example, for a < 0.015 we require a' < 0.01, which is much stronger 
than the corresponding post-Newtonian bound. 
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Future observations will be able to constrain n s to within 0.1 |]37| which would further 
improve the bounds on (a, a'), in the context of a given inflaton potential. Although far 
harder to measure, the tilt of the tensor perturbations, n g , gives a model-independent bound 
on a. 
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FIGURE CAPTIONS 

Fig. 1. The evolution of the amplitude of density perturbations after Hubble-crossing for a 
range of comoving scales, k*, as a function of the number A of e-folds from end of inflation, 



for the model described in Sect. VI 



Fig. 2. Classical trajectory in the space of fields (x, y) defined in Eq. ( 6.7|) . The solid line 
corresponds to parameters (ai = 10~ 3 ,a2 = 10~ 2 ), while the dashed line corresponds to 
(a x = 10" 2 ,a 2 = 10~ 2 ). 

Fig. 3a. The solid line shows the spectrum of density perturbations at the end of inflation, 
V^ e (k), as a function of the number A* of e-folds before the end of inflation, when the 
corresponding scale left the horizon, for parameters (a± = 10~ 3 , a 2 = 10~ 2 ). The dashed line 
corresponds to T , Q 1 (k) and the dotted line to Vq 3 {k). 
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Fig. 3b. Same as in Fig. 2a, but for parameters (ai = 10 2 ,a 2 = 10 2 ). 

Fig. 4. The tilt n(k) of the spectrum of density perturbations, as a function of the number AT* 
of e-folds before the end of inflation, when the corresponding scale left the horizon. The solid 
line corresponds to parameters (ai = 10~ 3 , a 2 = 10~ 2 ), while the dashed line corresponds to 
( ai = 10" 2 ,a 2 = 10~ 2 ). 

Fig. 5. The ratio R(k) of tensor (gravitational waves) to scalar (density) perturbations, as 
a function of the number N* of e-folds before the end of inflation, when the corresponding 
scale left the horizon. The solid line corresponds to parameters (a\ = 10~ 3 ,a 2 = 10~ 2 ), 
while the dashed line corresponds to (at = 10~ 2 ,a 2 = 10~ 2 ). 

Fig. 6. Contour plot for the tilt n s of the spectrum of density perturbations in the (a 1; a 2 ) 
parameter-space. The region below the curves is the allowed region for n > 0.6, 0.7, 0.8 and 
0.9, from top to bottom, (ai, a 2 ) corresponds to (a, a') at the end of inflation. For compar- 
ison, the dashed line corresponds to the post-Newtonian bounds on (a, a') in Eq. fl2.6Q . 

Fig. 7. Contour plot for the ratio R of tensor to scalar perturbations in the (a 1; a 2 ) parameter- 
space. The region below the curves is the allowed region for R < 4, 2, 1 and 0.5, from top 
to bottom, (ai, a 2 ) corresponds to (a, a') at the end of inflation. 
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